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A new charge quantization in a phase-polarized Cooper Pair Pump (CPP) is proposed, based 
on the topological properties of its Hamiltonian ground state over a three-dimensional parameter 
space P. The charge is quantized using a set of path in P covering the surface of a torus, and is 
a multiple of the integer Chern index ci of this surface. This quantization is asymptotic but the 
pumped charge converges rapidly to the quantized value with the increase in the path frequency. 
The topological nature of the current makes this CPP implementation an excellent candidate for a 
metrological current standard. 

PACS numbers: 85.25.Cp, 03.65.Vf, 74.50.+r, 74.78.Na 
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Topological defects occur in a number of different set- 
tings. In superfluids or magnets, the ground state is de- 
scribed by an order parameter. Topological defects (vor- 
tices, skyrmions...) are singular points of the order pa- 
rameter field [1]. A topological charge is assigned to the 
mapping from physical space to order parameter space. 
In quantum systems with few degrees of freedom, the en- 
ergy levels and eigenstates may depend on a set of con- 
tinuous classical parameters (or band indices). When the 
number of parameters is sufficient (three for a complex 
Hamiltonian), isolated degeneracies between two neigh- 
boring levels (or bands) can occur. The presence of a de- 
generacy has many physical consequences. For quantum 
wavefunctions, the degeneracies can then be viewed as a 
defect in their phase field [2]. Some classical examples are 
well known in polyatomic molecules [3], where the nuclear 
coordinates form a set of semiclassical parameters (within 
the Born-Oppenheimer approximation) for the electronic 
structure. Energy manifolds can have conical intersec- 
tions at isolated values of the nuclear coordinates [4, 5]. 
When a topological charge is assigned to these degen- 
eracies, the change of multiplicity of rotation-vibrations 
levels as a function of nuclear coordinates are easily un- 
derstood in topological terms. In molecular magnets, iso- 
lated degeneracies have also been found for specific direc- 
tions and values of the magnetic field[6, 7]. The integer 
quantum Hall effect has also been interpreted in term of 
Chern indices [8, 9] which are the total topological charge 
included within a filled Landau band. In all these ex- 
amples, a topological charge is associated to an isolated 
degeneracy between neighboring bands, where the defect 
in the wavefunctions phase field occur. Its physics can 
be described by a "monopole" placed at the degenerate 
point in parameter space. The flux of the magnetic field 
produced by the monopole through any surface enclos- 
ing the defect is proportional to the Chern index for this 
surface. For a the two dimensional electron gas, the Hall 



conductivity of a band is a xy = j x ci where c\ is the 
Chern index assigned to the band [8]. The phase change 
of a wavefunction on a closed path is proportional to the 
integral of the potential vector, which is also the flux 
through the surface bounded by the path. This phase 
change, better known as Berry's phase, takes a partic- 
ular value (tt) for some of the physically relevant paths 
enclosing degeneracies. In molecular magnets, this phase 
change leads to destructive interferences which quench 
quantum tunnelling at the degeneracy points. 

In this letter, the relevance of this physics to super- 
conducting circuits is illustrated around an example, the 
Cooper Pair Pump (CPP), which lowest two energy levels 
form two bands which depend on continuous parameters 
(gate charges or voltages and magnetic flux). The param- 
eter space is in this case three dimensional and degenera- 
cies (defects) occur at isolated points in parameter space. 
We demonstrate how a charge quantization, proportional 
to the defect topological charge, can be exploited to pro- 
duce a current source with metrological accuracy. 

The CPP-circuit represented in Fig. 1 includes a small 
inductance in serie with the Cooper pair pump. The CPP 
has three Josephson junctions in series which define two 
superconducting islands. As long as the device is domi- 



nated by the Coulomb charging energy Eq 
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islands, the topological properties of the circuit do not 
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FIG. 1: The Cooper Pair Pump. 
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FIG. 2: (Color) Plot of the energies (in units of the charg- 
ing energy) of the two lowest levels as a function of the gate 
charges n gs and n g d at ip = tt in the vicinity of the vertex Ti. 
This "diabolical point" is a conical intersection of the energy 
sheets between the two lowest bands. It represents a topo- 
logical defect in the wavefunction phase fields and cannot be 
removed by any perturbations which merely shift its position 
in parameter space. 



depend on the precise value of the Josephson couplings 
Ej or capacitances C, which can be taken as identical. 
Let rii and ri2 be the number of Cooper pairs in excess 
on each island and n g i = 9 \ e al •> n g 2 = 9 2 e g2 •> the po- 
larization charge induced by two gate voltages V g i and 
V g 2- The electrostatic charges on each island are n\ — n g i 
and ri2 — n g 2, where n g \ and n g 2 are classical variables 
tuning the islands electrostatic energies. The inductance 
L serves two purposes: the phase <$> across the CPP is 
controlled by the magnetic flux $ it threads, and it is 
also used as the input for the readout circuit measuring 
the circulating current through the pump. For a small 

inductance, the magnetic energy (^) — 2tt-^J defines 

a sufficiently deep potential minimum: quantum fluctu- 
ations of (j) across the CPP are quenched and the phase 
bias is set by the "classical" phase ip = 27rjk The circuit 
energy can be tuned through three independent parame- 
ters n g i , n g 2 and ip defining the parameter space P. 

Charging Tic and Josephson 7ij Hamiltonian con- 
tribute to the CPP energy. It is convenient to write the 
charging energy in term of h s = hi +fi2, the total charge 
on the CPP and hd = hi — fi2, the charge imbalance 
between the two islands. Their canonical conjugates, the 
phase variables O s and enter in the Josephson Hamil- 
tonian. In term of h s and hd the charging Hamiltonian 
is 



7~Lc = Ec 
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where E c = n gs = n gl + n g2 , n gd = n gl - 

are the corresponding polarization charges and — 3 
is a capacitance ratio. Each eigenstate \n s ,rid) is the 
ground state of He inside an hexagonal area ft(n s ,n d )> 
centered at the point (n gs = n Sl n g d= rid), in the n gs -n g d 
plane (shown as the base plane in Fig. 3). Its boundaries 
are electrostatic degeneracy lines between "neighboring" 
charge states, and the vertices are points of triple de- 
generacies. The vertices {Ti} = (n gs = ^n g d — 0) 
and {T2} = {rigs = \^ n gd — 0) form the unit cell of this 
hexagonal lattice. With an appropriate gauge choice, the 
Josephson Hamiltonian TLj depends on the phase bias ip 
as 

Hj = -2Ej cos 6 S cos G d - Ej cos(2G d + cp) . (2) 

The total Hamiltonian is 27r-periodic in ip. Lattice trans- 
lations in the hexagonal lattice (n gs -n g d) connect equiva- 
lent but physically different states. Because the parame- 
ter space P is three dimensional, TLj lifts the electrostatic 
degeneracies except at isolated points where the ground 
level is still degenerate, in P (a consequence of the von 
Neumann- Wigner theorem[10]). These points occur at 
the vertices {Ti} and {T2} of the hexagonal lattice for 
ip = 7r, as illustrated in Fig. 2. In realistic implementa- 
tion of the Cooper pair pump, it is favorable to keep the 
ratio between Josephson and charging energy f3 = 
small, an assumption made throughout this paper. 

In the limit of zero Josephson coupling, the electro- 
static states |0, 0), |1, 1) and |1, — 1) are degenerate at the 
vertex {Ti}. For finite Josephson couplings, the ground 
state remains degenerate at this point when cp — tt. The 
two degenerate ground states at Ti = {Ti,7r} are 
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with energy E g (T 1 ) = § E c 
excited state 
state {E e = 



Ej while the first 



oR = (a 



n r 



e) is well separated from the ground 
Eg + ^rp-). For a small deviation 
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n gd ^ 



tt) from degener- 
acy Ti, a two-level approximation (\g±)) is appropri- 
ate for small (3. The eigenenergies E±(SH) — E g db 

Y^P 2 ip 2 display a characteristic conical 
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intersection at Ti . After rescaling the deviation dR from 
Ti, as b x = § E c cr, b v = Ej tp,b z = ^ 5, the two- 
level Hamiltonian maps onto an isotropic spin-^ problem, 



'^±(T 1 )(R) 



b(R). 



(4) 
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This makes the conical nature of the intersection explicit 
(E±(SR) — E'pzb-^y-), and spinor eigenstates can be spec- 
ified using the direction of b rather than 5H. 

The topological nature of this degeneracy becomes 
clear when considering the adiabatic evolution of quan- 
tum states along different closed paths T "around" this 



3 



R (T) = R (0) 



> 





















6S (strip) 




FIG. 3: Representation of the degeneracy point {Ti, <^ = tt} in 
the parameter space P. The cylindrical surface which encloses 
this point has the topology of a torus. The Chern index is 
defined as the flux through this surface normalized by 2ty. 
Several paths on this cylinder are considered in the text. 



point. One of them is represented in Fig. 3 as the time 
evolution of a point R(£) in P. If \g±(H(t))) are instan- 
taneous eigenstates of the Hamiltonian, after a period 
T around any closed path \g±(T)) = e~ % ( r]± ~ lf± " ) \g±(0)) 
depends on two phases, rj± = \ J Q E g± (H(t))dt, the dy- 
namical phase and j± = i <f r A±(R) -dR, the geometrical 



phase. It is expressed as an integral of the vector poten- 
tial A±(R) = (#±(R)|V#±(R)) over the path. If the 
closed paths considered are sufficiently close to the de- 
generacy Ti where the two-level approximation is valid, 
Berry's geometrical phase can also be computed using the 
spin- 1 representation as j± = i <f F A±(b) • db, where the 
vector potential is azimuthal in this textbook example: 
A±(b) = ± 2 C |b S | sSfl e ^ Euler angles 0, (j) specify the 
direction of b). Its line integral over a closed path V is 
7±(r) = T^p-? where Sl(T) is the solid angle seen from 
the degeneracy. The vector potentials A± (b) can be rec- 
ognized as the gauge field of Dirac monopoles of strength 
=F^ placed at TV They produce a radial magnetic field 
B± = T27B13" which flux is ^2tt when integrated over any 
surface enclosing this single degeneracy. When normal- 
ized to 27r, this flux gives the Chern index c^ ±} = =pl 
of the surface with respect to the eigenstates \g±) and 
one assignes to the degeneracy Ti topological charges 

</:(T,) =1. 

This topological picture is not restricted to the im- 
mediate neighborhood of Ti where a two-level approx- 
imation can be made. Let us consider the cylindrical 
surface represented in Fig. 3. Since the top and bot- 
tom planes ip = and cp = 2tt are physically equivalent, 
this cylinder has the topology of a torus, and its surface 



is the constant radius |R| = p cylinder. The magnetic 
flux of B_ = V x A_ through this closed surface de- 
fines its Chern index c[~ ] of the ground state |#_(R))[11]. 
This index is non-zero only if the topology is non trivial, 
i.e. if the surface encloses degeneracies. Here the "cylin- 
drical surface" encloses only the isolated degeneracy Ti 
and Gauss theorem guarantees that its index is identical 
to the one of the small sphere around Ti just consid- 
ered. For a cylinder around T2 = {T2,7r} the Chern in- 
dices have opposite sign and q±(T2) = ±1. For a larger 
surface, the Chern index is the sum of the topological 
charges it encloses and when no degeneracies are present, 
it is zero. 

In a CPP, pumping of charges through the circuit is 
achieved when moving the ground state |g_(R)) adia- 
batically along a one-dimensional path in P. In the liter- 
ature, the paths considered so far have been closed tra- 
jectories around {Ti} or {T 2 } at constant cp for which 
the charge transferred is not quantized. When the path 
is chosen to cover a surface enclosing one of the singular- 
ities, the charge transferred becomes an integer multiple 
of the ground state Chern index. This quantization is not 
only observable but can be exploited to make a current 
standard of metrological accuracy, in a similar way as the 
quantum Hall effect is a standard of resistance based on 
the Chern indices of Landau bands. Let us consider a 
helical path r^ on the surface of the cylinder (Fig. 3) 
starting from the plane cp = cpo and ending in the plane 
P — ^0 + 27r while making an integer number p of wind- 
ings (in Fig. 3, ipo = and p = 2). The pumped charge 
Qpump on pg^h a dynamical Q dyn and a geomet- 
rical contribution Q geo . The expressions 



ldyn 



), Q S eo = _ 2e ^ (r )f (5) 
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proved elsewere[12], generalize earlier results[13, 14] to 
arbitrary three dimensional closed paths in P. We first 
show that the geometrical contribution is quantized when 
averaged over the initial phase cpo of the helix. Taking 
two helices and r^ 0+( ^ shifted in (p by an infinites- 
imal increment Sep, a closed path S(y?o) on the surface 
of the cylinder (shown in Fig. 3) can be constructed by 
adding two small opposite "vertical" segments connect- 
ing the paths r^ and r~* +5 which line integrals cancel. 
Using Eq. 5, the transferred charge is 
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geo 



2e 



7-OVo) 



S(p 



o<p 



Using Stokes theorem, the above numerator is also 
J 55 B_(R) • hdS, where SS is the surface of the strip 
between the two helices (shown in Fig. 3 for tpo = 0) and 
B_(R) is the "magnetic" field produced by the monopole 
at the degeneracy. When the initial angle tpo is integrated 
from to 27r, adding all the elementary strips 5S gener- 
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FIG. 4: Plot of the peak amplitude of Q dyn (r v?0 ) as a function 
of the number of windings p in the interval [(fo ; ipo + 2tt] for a 
charging energy of lOmeV, a Josephson coupling Ej = 0.5 Ec 
and a period of 10 -6 sec. 

ates p times the total surface S (Fig. 3). Hence 
(Q Se °) V o = § (£b_ • hdS = 2epc[-\S) =2ep (7) 

is quantized by the Chern index c\ of the cylindrical sur- 
face S (torus). It is not necessary to average over ipo to 
reach quantization. If the number of windings p in T is 
large, the path integral of A_ over T reaches the sur- 
face integral of B_ over the cylinder and Q geo p ^=°° 2ep. 
The increase of accuracy with p being roughly an or- 
der of magnitude per additional winding, this is not only 
practical but also quite accurate. Using the same proce- 
dures, the dynamical charge Q dyn averages out to zero 
because 7^_(r^ ) is a 27r-periodic function of cpo since the 

Hamiltonian is 27r-periodic in cp. Its magnitude ~ 

depends critically on the ratio /3 = J^, however the pe- 
riodic and odd-dependence of Qdyn with cpo guarantees 
that its (/?o-average is zero. In realistic situations, (po may 
not be completely controlled, or the path may not close 
perfectly because of flux- noise. It is therefore instructive 
to plot the evolution of the peak amplitude of Q dyn as a 
function of the number of winding p in a ip— period. In 
Fig. 4, this peak amplitude looses four order of magni- 
tudes when p is increased by 10, using a ratio (3 = 0.5 
an order of magnitude larger than its optimal value. For 
smaller /?, the decrease is even more dramatic. Hence, 
the CPP delivers a very accurate DC current 

jpump gpump y V^oo ^ ^ ^ ^ 

where v is the winding frequency. 

We now consider the possible sources of errors in a 
practical implementation of this circuit as a current stan- 
dard. Low frequency noise in the parameters will intro- 
duce some jitter in the path covered in P: the surface gen- 
erated by the path will not be a perfect cylinder. On the 
other hand, as long as the jitter amplitude is small com- 
pared to the radius, the surface generated still encloses 



the topological charge, and no significant error in quanti- 
zation should result. The most detrimental source of er- 
ror comes from inelastic transitions between the ground 
state and the first excited state in the adiabatic evolu- 
tion, because their topological charges are opposite. This 
means that if the time spent in the first exited state is 
r in a period, the relative error in the transferred charge 
will be — 4e^r. There are two ways this can happen. For 
each winding around T\ , three saddle points between the 
ground and first excited manifolds have to be crossed 
(along the hexagonal lines shown in Fig. 3), where the 
probability of Landau- Zener tunnelling is largest. For 
the optimal radius (p± |) this tunnelling probability is of 

/ 3?r \2 Ej Ej 

order P ~ e ^ 2 ' E c h » . To stay below the part per mil- 
lion level requires that the winding frequency v to stay 
below 100 MHz for a Josephson coupling of 0.05i^cs lim- 
iting the pumped current amplitude 2ev to 32 pA. The 
other source of inelastic transitions to the first excited 
state is the back-action of the measuring device on the 
CPP, which needs to be rigourously controlled. 

The physics discussed in this paper is by no means lim- 
ited to the CPP: as long as a circuit has three tunable 
parameters or more, degeneracies can occur. In this pa- 
per, we emphasized the topological quantization and its 
strong robustness to adiabatic parameters fluctuations, a 
key point for metrological applications. The realization 
of this metrological source is under way in our group: it is 
a beautiful scientific challenge for fundamental quantum 
electronics of real practical value. R. Leone is supported 
by a fellowship of the Rhone- Alpes Micro-Nano cluster. 
The authors are grateful to F. Faure for enlightening dis- 
cussions. 



[I] N. D. Mermin, Rev. Mod. Phys. 51, 591 (1979). 

[2] M. V. Berry, Proc. R. Soc. Lond. A 392, 45 (1984). 

[3] G. Herzberg and H. C. Longuet-Higgins, Disc. Faraday 

Soc. 35, 77 (1963). 
[4] F. Faure and B. Zhilinskii, Phys. Rev. Lett. 85, 960 (2000). 
[5] B. I. Zhilinskii, Phys. Rep. 341, 85 (2001). 
[6] W. Wernsdorfer and R. Sessoli, Science, 284, 133 (1999); 

W. Wernsdorfer, N. E. Chakov and G. Christou, Phys. 

Rev. Lett. 95, 037203 (2005). 
[7] P. Bruno, Phys. Rev. Lett. 96, 117208 (2006). 
[8] D. J. Thouless, M. Kohmoto, P. Nightingale and M. den 

Nijs, Phys. Rev. Lett. 49, 405 (1982); D. J. Thouless, 

Phys. Rev. B 27, 6083 (1983). 
[9] M. Kohmoto, Ann. of Phys. 160, 296 (1984). 
[10] J. von Neumann and E. P. Wigner, Phys. Zschr. 30, 467 

(1929). 

[II] B. Simon, Phys. Rev. Lett. 51, 2167 (1983). 
[12] R. Leone and L. P. Levy, (unpublished). 

[13] M. Aunola and J. J. Toppari, Phys. Rev. B 68, 020502(R) 
(2003). 

[14] M. Mottonen, J. P. Pekola, J. J. Vartinainen, V. Brosco, 
F. Hekking, Phys. Rev. B 73, 214523 (2006). 



